Propagation of non-classical correlations across a quantum spin chain 
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We study the transport of quantum correlations across a chain of interacting spin- 1/2 particles. As a quanti- 
tative figure of merit, we choose a symmetric version of quantum discord and compare it with the transported 
entanglement, addressing various operating regimes of the spin medium. Discord turns out to be better trans- 
ported for a wide range of working points and initial conditions of the system. We relate this behavior to the 
efficiency of propagation of a single excitation across the spin chain. Moreover, we point out the role played 
by a magnetic field in the dynamics of discord in the effective channel embodied by the chain. Our analysis 
can be interestingly extended to transport processes in more complex networks and the study of non-classical 
correlations under general quantum channels. 

PACS numbers: 03.65.Yz, 75.10.Pq, 42.50.Lc 
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The behavior of features such as quantum coherence and 
entanglement in a composite quantum system whose state is 
exposed to the effects of environmental actions has been the 
focus of an extensive research activity. Recently, much atten- 
tion has been paid to the case of environments embodied by 
systems of interacting quantum particles [1]. Such dynami- 
cal environments can induce interesting back actions on the 
evolution of a system, thus significantly affecting its proper- 
ties. From the point of view of coherent information process- 
ing, on the other hand, the non-trivial dispersion properties of 
networks of such interacting particles represent an interesting 
opportunity for their use as short-haul communication chan- 
nels for the inter-connections among on-chip nodes in the next 
generation of information processing devices H. 

While most of the work in these contexts has focused on 
the study of the properties of entanglement upon propagation 
in such media, it is now widely accepted that the space of non- 
classical correlations accommodates more than just quantum 
entanglement. Figures of merit such as quantum discord [3, 4] 
and measurement-induced disturbance D5|], to cite only two of 
the most popular ones, are able to capture the content of non- 
classical correlations of a state well beyond entanglement. Al- 
though the role played by such broader forms of non-classical 
correlations in the quantum mechanical manipulation of infor- 
mation has yet to be fully understood, enormous is the interest 
they bring about as the manifestation of the various facets of 
quantumness in a system. It is thus very important to work 
on the exploration of the behavior of such quantities upon ex- 
posure to dynamical and finite environments of the sort ad- 
dressed above, so as to build a useful parallel with the much 
more extensively investigated case of entanglement. 

In this paper we study the propagation of quantum corre- 
lations across a system of interacting spin-1 /2 particles. Our 



main goal is to compare the way important indicators of non- 
classicality, such as quantum discord (QD) JH 0] and en- 
tanglement of formation (EoF) |6|. are transferred through a 
medium offering non-trivial dispersion properties. In doing 
this, we aim at understanding whether or not the fundamen- 
tally conceptual difference between entanglement and discord 
leaves signatures in the way such non-classical quantities are 
transferred. We show that this is indeed the case by preparing 
a non-separable (in general mixed) state of an isolated spin 
and the one occupying the first site of a linear spin-chain. We 
then compare the quantum-correlation properties of such an 
initial state with those of the state achieved, at a given instant 
of time of the evolution, between the isolated spin and the one 
occupying the last site of the chain itself. QD appears to be 
better transmitted than entanglement (as quantified by EoF) 
in a wide range of working conditions and regardless of the 
details of the initial state being considered. It is more robust 
to the dispersion inherent in the effective spin-medium across 
which it propagates, being non-zero in situations where the 
EoF is, for all practical purposes, null. By relating the en- 
tanglement to the single-excitation transition amplitude of the 
system, we identify the working point at which a cross-over 
occurs between the quality of transport of QD and EoF, mak- 
ing the transport of entanglement more efficient. Moreover, 
interesting effects of entanglement-forerunning [7], where QD 
precedes the establishment of EoF, are found in the way quan- 
tum correlations build up between the isolated spin and the 
last one in a given chain. Our analysis considers a large num- 
ber of state families in such a transport problem, addressing 
explicitly those that maximize the degree of discord at given 
global mixedness ifsj- TToll . 

Our study provides exact quantitative answers to a problem 
that has been so far largely overlooked, although being rele- 



vant for a wide range of physical situations. For instance, it 
is sufficient to think about recent studies of the propagation 
of information in biological systems operating on the verge 
of quantumness, which a ppe ar to benefit from the inclusion 
of a mild degree of noise I! ill . In such conditions, one might 
wonder whether other forms of quantum correlations are fa- 
vored, given that the perfect transport of entanglement would 
be prevented. 

The remainder of this paper is organized as follows. In 
Sec. |T]we describe the physical situation at hand and provide 
the general analytical form of the time-dependent density ma- 
trix describing the state of the isolated spin and the last one 
in a chain. Sec. [TT] is devoted to a brief introduction to the 
quantitative indicators of quantum correlations adopted in this 
work. In Sec. IIII Al we face the propagation of EoF and QD 
across the spin chain and address the relation between such 
non-classicality indicators and the single-excitation transition 
amplitude. While the transport of discord appears to be fa- 
vored, we point out the existence of a cross-over point in the 
parameter space at which the performance of EoF becomes 
superior. Our study includes both pure and mixed input states, 
among them the case embodied by states that maximize QD 
at fixed global entropy. Sec. IIII B I addresses the effect of a 
uniform magnetic field on the dynamics of the discord find- 
ing that it has relevance only for states lacking of rotational 
symmetry in the XT-spin plane and, in Sec. IIII CI we describe 
briefly the evolution of the non-classicality indicators under 
the influence of an environment modelled by a spin chain. Fi- 
nally, in Sec. [IV] our conclusions are drawn and some open 
questions arising from the present work are put forward. 



I. THE MODEL 

We consider the configuration shown in Fig. \T\ i.e., a quan- 
tum channel consisting of N interacting spin- 1/2 particles in 
a linear configuration with open boundary conditions. The 
Hamiltonian model describing the system is taken to be (we 
take units such that ft = 1 throughout the manuscript) 



H 



<fi = -27 Y^S'J'; 1 +S i y S i y +l )-2hY j S i z 



(1) 



with S' k the k-x,y, z spin-component operator of particle 
i = 1, ...,N, J the inter-spin coupling strength and h a uniform 
magnetic field. In what follows, the spins occupying sites 
j = 2, ...,N will be assumed to be all prepared in down state 
IX), with {|X) , |t)} denoting the eigen-states of S z . On the other 
hand, the first spin forms a (generally mixed) bipartite quan- 
tum correlated state with a further particle, labelled 0, which 
is physically detached from the chain. 

The symmetry properties of the Hamiltonian model in Eq. Q] 
restrict the dynamics to those states within the zero- and 
the single-excitation sectors of the total Hilbert space of the 
spins 111 2n, where the chain behaves as an amplitude damping 
channel lfl3ll fully characterized by the transition amplitude of 
the spin excitation from site 1 to site r= 1, ...,N. The latter 
is conveniently expressed as f r (t) = {r\ e~ fH ' |1), where in the 
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FIG. 1: (Color online) We consider a chain of A' interacting spin- 1/2 
particles coupled through the Hamiltonian model "H [cf. Eq. (J)]. A 
further spin- 1/2 particle, labeled and completely isolated from the 
chain, is prepared in a joint mixed state p (10) (0) with particle 1. We 
study how the general quantum correlations of such a state propagate 
across the chain. 



states \n) (n = 1, N) all of the spins are in |X) except the one 
at position «, which is in |f). 

In this basis, *K is represented by an N x N tridiagonal ma- 
trix, which can be analytically diagonalized for any length. 
The reduced density matrix p (r) (f) describing the state of par- 
ticle r at time t can thus be expressed as a function of the state 
of particle 1 at the reference time t — 0. That is 



p^(0)+p^(0)(l 



\fr(t)\ 2 ) 



p;- (o)/;w 



P$(0)|/ r (*)l 2 



, (2) 



where p a0 = (a\p\p) (a,/3= T,T)- 

In the remainder of this work we will consider the chain as 
prepared in the fully factorized state |X)j, while the initial 
state of spins and 1 is a quantum correlated state. Clearly, 
while spin undergoes only a free evolution, the chain's ele- 
ments evolve according to the intra-chain coupling model in 
Eq. (Q~|). This implies that the overall time evolution operator 
*ZY(f) that propagates the state of the N + 1 spins factorizes as 
f U(t) — i ® e~'' Ht . Together with the uncorrected initial state 
of the remaining part of the chain (which thus shares no corre- 
lation with spins and 1), this legitimately allows us to make 
use of the formalism developed in Ref. IU411 to get the joint 
state p (r0) (f) of spins r and from the knowledge of p (r) (f)- 
Notice also that a very similar approach, which has also been 
verified by an exact numerical study, has been previously used 
in similar contexts lfl5ll . Although such an approach is broadly 
valid and can indeed be used for any initial state of spins 
and 1, here we restrict our analysis to X-type input states of 
the general form 
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(3) 
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Here we have introduced the compact notation 
Pij=(i\p (m (0) \j) with |l) = |XX) , \2) = IXT) , |3> = IU> , |4> = 
ITT)- in f act ' as mentioned above, the excitation-preserving 
nature of the Hamiltonian studied here (which commutes 
with the total number of excitations in the system) ensures 
that the X-type character of any input state of the pair (1 , 0) is 
preserved upon evolution. More explicitly, the only non-zero 



3 



elements of the evolved density matrix will be 



p\ r f\t) = PlI + (1 - |/,.(r)l z )P33 , p£ u '(r) = I/XOIW 

P^f(f)=P22 + (l-|/ f -(f)| 2 )P44 , P^ 0) (0 = l/r(0| 2 p44, (4) 
P ( lf (t)=fr(t)PH , P^ 0) (0 = / r (f)P23. 

The initial conditions being specified by the input state, the 
above equations describe both the propagation of quantum 
correlations from site 1 to site r, when r > 1 , and the decoher- 
ing influence of an environment (embodied by the spin chain) 
on the spin occupying site 1, when r — 1. In the former case, 
the single-excitation transition amplitudes f r (t) contain infor- 
mation on the working conditions of the channel. The limiting 
case f r (t) = 1 gives p < - r '°\t) =p (1 ' 0) (0), i.e., the perfect transfer 
of the input state from pair (1 , 0) to (r, 0). On the other hand, 
fx (f) defines the probability amplitude of finding the excitation 
on the first spin. Note that X-type density matrices are such 
that no single-spin coherence will develop in time, i.e., both 
the reduced single-spin density matrices remain diagonal. 
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II. FIGURES OF MERIT FOR QUANTUM 
CORRELATIONS 



After having introduced the dynamical model that will be 
addressed in our study, we turn our attention to the figures 
of merit that will be used in order to perform our quantita- 
tive analysis. As already anticipated, we take QD yllj] as a 
measure for general quantum correlations between any two 
spins under study. As originally proposed by Ollivier and 
Zurek, QD can be associated with the difference between two 
classically equivalent versions of mutual information, which 
measures the total correlations within a quantum state. For a 
two-spin state p (r ' r ' extracted from our system, the mutual in- 
formation is defined as J(p (r - r ' ) ) = S(p (r) )+S(p {r ' ) )-S(p {r ' r ' ) ). 
Here, S(p) = -Tr[p log 2 p] is the von Neumann entropy of 
a generic state p. Alternatively, one can consider the one- 
way classical correlation J^(p (r ' r ' ) ) = S(p (r) )-'H [ f ll] (r\r') fl, 
where we have introduced < H^(r\r')= £j piS(p' r ^,) as the 
quantum conditional entropy associated with the the post- 
measurement density matrix p'^ r , = TT r >[flip^ ,r ) ]/pi obtained 
upon performing the complete projective measurement {II,} 
on spin r' (p/ = TrflTp^ ']). QD is thus defined as 



)-J* 



(5) 



with the infimum calculated over the set of projectors {ft,} 
(3D. Analogously, one can define £)""*, which is obtained upon 
swapping the roles of r and r' . The inherently asymmetric def- 
inition of QD makes, quite naturally, +Tf~ '. This might 
be the cause of misinterpretations: a quantum-classical state 
for which Tt~ + but !D~~* = (or viceversa) lfl6ll might be 
interpreted as strictly classical if only (!D*~) is probed. 
Here, we are interested in the transport of quantum correla- 
tions, regardless of the way they are encoded in the two-spin 
state. Rather refined solutions to this issue passing through 
the generalization of the definition of QD or the introduc- 
tion of strictly faithful entropic measures that are null only for 



classical-classical states (i.e., states such that £) <_ ' _> =0) have 
been proposed [9]. However, they typically require a double 
optimization to be performed over a bilateral set of projective 
measurements. In order to bypass the numerical burden that 
this would imply, we consider the two-way QD, 



£> = max[£> < ", £>"*], 



(6) 



which is strictly null only on states endowed with no quan- 
tum correlations, and faithfully signals classical-classical 
states fH]. 

On the other hand, our chosen entanglement measure is 
EoF (6(], which quantifies the minimum number of Bell pairs 
needed in order to prepare a copy of the state p ( ' v ' ' we are 
studying. The relationship between EoF and QD has been 
recently examined to study the distribution of quantum corre- 
lated states in the entropic space BSD- Moreover, it is possi- 
ble to establish a triangular relation connecting QD, EoF and 
conditional entropy in multi-spin quantum states [ 17], so that 
such two figures of merit appear to be natural choices for a 
quantitative comparison. For arbitrary two-spin states, EoF is 
calculated as 



£ = /z^[l + Vl - C 2 ] 



(7) 



where h(x) = - xlog 2 x - (1 - jc)log 2 (l - x) is the binary en- 
tropy function and C is the concurrence of the state |@]. The 
latter, an equally valid entanglement measure, is found in 
terms of the eigenvalues A\ >/l2,3,4 of the matrix p^ ,r '\fr y <g> 

&y)p < - r,r '^*(&y ® fry) 8S 



C = max 



4 

o, 



(8) 



where fr y is the y-Pauli operator. For an X-type state of pair 
(r, r'), the concurrence is straightforwardly shown to be 

C inr r } = 2 max [0, |p 14 | - VP22P33, IP23I - VP11P44J ■ (9) 

Differently from EoF, QD does not have a closed analytical 
expression for any two-spin state, although some steps toward 
this goal have been performed |[l8l[l9ll . Nevertheless, for the 
special case of the class of states presented in Eq. ©, one can 
obtain analytic formulas for D as a function of the dynamical 
parameter /(f), as well as of the input state. However, as their 
expressions are lengthly and not very informative, we do not 
report them explicitly here. 



IE. PROPAGATION OF QUANTUM CORRELATIONS 

A. Case Study of Pure States and MDMS 

We perform our analysis by addressing the propagation of 
quantum correlations across the spin chain when the pair of 
spins (1,0) is initialized in a given non-separable state, while 
the remaining spins are in ®^ =2 In order to address the 
temporal behavior of our figures of merit, we need the explicit 
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form taken by the single-excitation transition amplitude f r (t). 
Although the present formalism is valid without any major 
difficulties for general r belonging to the chain, for the sake 
of clarity we set hereafter r-N and omit the subscript in f(f). 
For a uniform chain ruled by Eq. (Q3 such quantity reads 



(10) 



k=\ 



We will use this explicit result to build up the state of spins 
and N and thus calculate quantum correlations. The single- 
excitation transition amplitude is a real (purely imaginary) 
quantity for N odd (even), due to the symmetry properties of 
the spectrum of the system (a more detailed discussion on the 
properties of function /(f) is given in Ref. J28ll ). 

As a first significant instance, we consider the case in which 
the joint state of spins and 1 is pure. In fact, for pure states 
£ = T> and it is interesting to study whether or not, in this 
case, entanglement is lost in favor of discord as information 
propagates across our dispersive medium. As the EoF is based 
on concurrence, we consider the class of pure entangled input 
states parameterized as 

k(Qi,0))> (1 Q) = siny |U> + cos y |TT> , d D 

where y=(l/2) arcsinC^o)- The evolved state is then achieved 
using the approach outlined in Sec. Q] 

Perfect transfer of entanglement across an interacting spin 
chain is known to occur under proper conditions I120ll . In 
particular, it was proved that perfect entanglement transfer is 
achieved when perfect end-to-end state transfer is made pos- 
sible. For a system of three spins, the Hamiltonian model in 
Eq. [T]allows for perfect state transfer 1I2TI1 . thus implying that 
this is also the case for entanglement and QD, since for pure 
states £ = £). Nonetheless, the examination of Fig. [2] shows 
some interesting features. Clearly, the peaks shown in all pan- 
els of Fig. |2] are achieved at the instants of time at which the 
input state dTTb is perfectly transmitted, where the pair (3 , 0) 
is pure and £(i,o> However, between two consecutive 
peaks, the state of spins and N is mixed and the two fig- 
ures of merit can be quantitatively different. Evidently, within 
these time windows the transport of QD is favored with re- 
spect to EoF, being not only quantitatively larger than £ but 
also non-null at times such that the EoF is, for all practical 
purposes, zero. Interestingly, as we increase the initial degree 
of entanglement of I^C^o))),, 0) , this effect becomes less im- 
portant until, at C(i_o) — 1, there are narrow regions close to 
the peaks where the transported entanglement overcomes QD. 
We will see later on in this paper how this effect depends on 
the single-excitation transition amplitude . 

The behavior disclosed above becomes much more visible 
when we consider longer chains. The uniformity of the cou- 
pling strengths across the medium prevents perfect state trans- 
fer (and thus perfect entanglement transfer) for chains of more 
than three spins 1I20I I2TI1 . Therefore, we find an interesting 
feature in the corresponding propagation of T> [see Fig. 0: 
the discrepancy between the propagated EoF and QD be- 
comes quantitatively much more significant, while the degree 
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FIG. 2: (Color online) Time behavior of quantum correlations be- 
tween spins and 3 in a system ruled by "H with uniform interac- 
tion strengths and arbitrary h, The panels (a)-(d) are for C(i,o) = 
0.1, 0.4, 0.7, 1 respectively, and the initial EoF is shown as a straight 
line. The solid (red) curve shows the EoF, while the dashed (gray) 
one is for the shared two-way QD. 



of propagated discord is damped in time much more slowly 
than entanglement. However, such effects depend strongly on 
the initial value of entanglement in a way that the behavior of 
£ and D almost merge as £(i,o) — > 1. In Fig. [3] we show an 
instance of this case by reporting the propagation of £ and T> 
across a system of 50 spins interacting according to 

Needless to say, as the two figures of merit refer to two dif- 
ferent forms of quantum correlations, some quantitative dif- 
ferences should be expected. However, here we would like to 
stress that, the input state pure, as remarked above, QD and 
EoF are exactly equivalent. This implies that the differences 
\D(t) - £(0)| and |£(f) - £(0)| faithfully quantify the perfor- 
mance of each non-classicality indicator upon propagation. 

These issues are better discussed by looking at Fig.|4] where 
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FIG. 3: (Color online) Time behavior of quantum correlations be- 
tween spin and 50 for a system ruled by "H T . with arbitrary h. The 
panels (a)-(d) are for Qi.o) = 0.1, 0.4, 0.7, 1 respectively, and the ini- 
tial EoF is shown as a straight line. The solid (red) curve shows the 
EoF £, while the dashed (gray) one is for the shared two-way QD D. 
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FIG. 4: (Color online) Behavior of QD against entanglement and propagation time for a chain of 15 spins, homogeneous intra-chain couplings 
and arbitrary h. The yellow plane at D=& is used as a guide to the eye for discerning whether or not D>&. Each panel is for a different 
initial value of concurrence in pair (1,0). We have taken C(i,oj = 0.1, 0.6 and 1 in panel (a), (b) and (c), respectively. As the initial degree of 
entanglement grows, the efficiency of transport of D gets very close to that for £. 



we plot the propagated QD and EoF against time for a total of 
15 spins. In each panel, the plane for which £(f) = £)(?), Vf 
is displayed as a reference. Any point lying above the plane 
corresponds to discord overcoming EoF. Our choice of N = 15 
is only due to the clarity of the associated figures, and does not 
hinder the validity of our conclusions. 

We now show that the qualitative features revealed by our 
study on pure input states hold also when a mixed state of the 
spin-pair (1,0) is prepared. As in general the value of dis- 
cord and entanglement associated to such input states will not 
coincide, this situation encompasses from the start the pro- 
found differences between QD and EoF. In order to wash out 
the ambiguities associated with the possible choices for input 
mixed states, we refer to the studies in Refs. H Ht]. There, 
families of two-spin mixed states maximizing the degree of 
the two-way QD at fixed values of the global entropy have 
been identified and fully characterized. Furthermore, for these 
states D*~ - TF* so it is immaterial which spin is attached to 
the chain. We dub such states as maximally discorded mixed 
states (MDMS), which represent the counterpart, as far as dis- 
cord is concerned, of the well-known maximally entangled 
mixed states (MEMS) f22ll . Even more interesting, under 
proper choices of entanglement measures, part of the MEMS 
frontier is shared with the MDMS one H. As such extremal 
states are clearly dressed with a particular significance, we 
now restrict our study to them. The MDMS boundary is a 




FIG. 5: (Color online) Comparison between the re-scaled quantities 
T) (grey dashed line) and £ (red solid one) propagated across our 
chain of uniform couplings and arbitrary magnetic field. The spin- 
pair (1, 0) is prepared in a Werner state with a = 0.4, 0.7 and 1 [panel 
(a), (b) and (c) respectively]. The curves are plotted against the tran- 
sition amplitude |/|. 



patch-work of three classes of X-type states whose elements, 
following the notation used in Eq. ©, are given by 
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a 




= Pl4 = 


2' 




1 + a 
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1 - a 




P R n 
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P22 = 



1 



a-g 



rw w 
P22 = P33 



2 

1 — a 



p 

p 33 = 



w 

Pl4 



1 — a + j 



(12) 



P22 



P?3 



Pi4- 2 



and p 23 ' = 0. While p is a general two-parameter family 
(a + g < 1), p w (-1/3 <a < 1) are Werner states JH and p R 
(0 < a < 1 /3 with a + g < 1) are MEMS when the relative en- 
tropy is chosen as a measure of entanglement [22]. Such states 
belong to the MDMS frontier only under properly chosen val- 
ues of a and g. Such conditions are in general highly non- 
trivial and we refer to Ref. [9] for full details. Here it is enough 
to state that we will consider values of such parameters that 
guarantee the MDMS nature of the corresponding states. 

In order to provide a faithful evaluation of their perfor- 
mance upon propagation, we will compare the re-scaled quan- 
tities fc = K/K(i,o) with = [%, )=£(i,o),£>(i,o)] being 
the value of one of our figures of merit after propagation [for 
the initial spin pair]. Moreover, rather than replicating the 
time-dependent study performed so far and in order to provide 
a universal analysis freed from the choice of N, we will con- 
sider the propagated QD and EoF as general functions of the 
single-excitation transition amplitude |/| e [0, 1] (from now 
on we drop the label stating its dependence on time). 

We start by studying Werner states, which are entangled 
only for a > 1/3. For values mildly larger than this threshold, 
where the purity of the state is small and also its entangle- 
ment, very large values of |/| are required in order to actually 
transport S. Differently, T> is non-null for any |/| and irre- 
spectively of the initial QD properties of pair (1,0). The rel- 
ative discrepancy between the two figures of merit is in gen- 
eral very large and decreases only for almost ideal transport 
of excitations across the chain. As a — * 1, i.e., by increasing 
the purity of the state, more EoF is transported, even in the 
low-/ region, thus reducing the differences between the two 
non-classicality indicators. In the limit of a = 1, which makes 
p w a maximally entangled pure state, discord is overtaken by 
the EoF at \f\ > 1/ V2. In fact, the state of pair (N, 0) corre- 
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sponding to such a value of the transition amplitude reads 



l 

2 





1 

2V2 





\ 



~~i= 



1 

2V2 



J_ 

4 



(13) 



which is an instance of a mixed state (its von Neumann en- 
tropy is 0.811278) with Fig. [5] exemplifies the cases 
discussed above for three different values of parameter a. 

When addressing the case of the two-parameter family p p , 
the situation is even more striking. This class of states is 
MDMS for g = and a € [0.503, 0.574], for instance, which 
correspond to a region of large mixedness such that the re- 
scaled QD is always larger than the corresponding re-scaled 
EoF, as seen in Fig. [6] In this case, S never overcomes T) and 
can only equal it at \ f\ = 1. 

Finally, the general picture is confirmed by the investiga- 
tion on the third family of boundary states, p R , which are also 
the states maximizing the relative entropy of entanglement 
at fixed global mixedness |22fl. The conditions that a and g 
should satisfy in order for p R to be a MDMS are rather non- 
trivial, passing through the solution of a transcendental equa- 
tion H. For our purposes, it is enough to state that p R spans 
the large-purity region of the boundary and thus the crossing 
of QD and EoF at some value of |/| should be expected. This 
is indeed the case, as highlighted in Fig. [7] for two values of 
parameters a and g that guarantee the MDMS nature of such 
class of states. While Fig. [7] (a) refers to a large-purity case 
(the von Neumann entropy of the corresponding state is 0. 159, 
implying a small degree of mixedness) where EoF wins over 
QD at large enough |/|, panel (b) is for a much more mixed 
state. In this case, as expected from the analysis above, T> is 
better transported at any dynamical condition. 



B. Effects of a magnetic field 

Here we study the effects of a uniform magnetic field h on 
the transport of quantum correlations. From Eq. (TTOi we see 
that the single-excitation transition amplitude at h + differs 




FIG. 6: (Color online) Comparison between the re-scaled quantities 
T) (grey dashed line) and £ (red solid one) propagated across our 
chain of uniform couplings and arbitrary magnetic field. The spin- 
pair (1,0) is prepared in state p^ |0) with b = and two values of a. 
The curves are plotted against the transition amplitude |/|. 




FIG. 7: (Color online) Comparison between the rescaled quantities 
T) (grey dashed line) and £ (red solid one) propagated across our 
chain of uniform couplings and no magnetic field. The spin-pair 
(1,0) is prepared in state p* with (a, g) = (0.0150, 0.9747) [panel 
(a)] and (a,g) = (0.1625, 0.7649) [panel (b)]. The curves are plotted 
against the transition amplitude |/|. 



from that at h — only by an overall oscillating phase factor, 
f <h \t) = e~' 2h 'f(°\t). Consequently, by using the general form 
of an evolved state in Eq. (0 and the expression for the con- 
currence of X-type states given in Eq. (0, it is straightforward 
to see that & would depend on just |/(f)| so that the introduc- 
tion of a magnetic field does not affect the propagation of EoF. 

The situation is radically different when considering the 
transport of QD for a certain class of input states of pair 
(1,0). In fact, by considering states with maximally mixed 
marginals H p (1 ' 0) (0) = \ (2 + E;= w c t &i ® <x;) (c, e R), for 
which analytic expressions for QD hold, we find that the mag- 
netic field increases the amount of discord that can be obtained 
between and N as compared to the case with h = 0. The rea- 
son for this enhancement can be found in the fact that, when 
\c x \ + \c y \, D depends on both the real and imaginary part of /. 
Moreover, the phase factor e~ aht in / ( ' !) (f) yields an oscillating 
behavior of D that is lower-bounded by the value achieved at 
h = 0. In Fig.[8]we show the typical behavior described above. 

X-type states with maximally mixed marginals allow for 
the identification of cases where the conditions for a truth- 
ful transport are breached. In all the cases studied so far, QD 
is a monotonically increasing function of \f \ with a maximum 
occurring at f* ~ N, which are all features consistent with 
the picture of a transmission mechanism. For a maximally- 
mixed-marginal state having \c x \ + cJ, on the other hand, the 
discord between spin and N can be larger than the initial 
QD content of p (I 0) . This point is best illustrated with the aid 
of an example. The density matrix elements of a two-spin X- 
type state with zero discord should satisfy one of the following 
conditions HIH 



1. P14 = P23 = 0. In this case, all the coherences are identically 
null and p is purely diagonal. The corresponding state is 
thus a classical-classical one as 2) <_ '^ = 0. 

2. pn = p22,P33 = P44and|pi4| = |p23|. This case corresponds 
to a quantum-classical state with !D~ > + and JT~ = 0. 

3. pn = P33, P22 = P44 and |pi 4 | = \p23l which give rise to a 
classical-quantum state with D - * = and !D*~ + 0. 

When the symmetrized discord in Eq. (O is used as a figure 
of merit, conditions 2 and 3 collapse into 
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4. p a =l/4 (z=l, ...,4) with |p 14 | = ^ 23 |. 

We now consider a p (I0) (0) having c y =c z = and c x =l, 
which corresponds to a state endowed with only classical cor- 



n\f\ 



relations (embodied by (fr x Qfr^') + 0). The evolution yields 
a quantum-classical state 
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If / is a real (imaginary) function, such state develops a non- 
zero correlation function ®<t^) ((fry ® o"x°')) that is 
quantitatively equal to / itself. However, spin N also de- 
velops a non-zero magnetization in the z-direction given by 
<<5-f° ) = 1 -|/| 2 . According to the results in Refs. dHHH, due 



to the non-commutativity between &f® ® i (0) and fr^ t 



-(()) 



the initially classical-type correlations acquire a quantum na- 
ture responsible for non-zero discord, as is shown in Fig. |9] 
Reasoning along the same lines, we can see that an initially 
zero-QD state fulfilling condition 1 remains such because all 
one- and two-spin correlators involve only fr z . 




22 24 26 28 30 32 34 



22 24 26 28 30 32 34 



FIG. 8: (Color online) QD transported across a chain of 50 spins 
whose (1,0) pair is initialized in an X-type state with maximally 
mixed marginals having c^-0.53, c y — 0.340, c,^0.035. The associ- 
ated value of QD is £> m (0) = 0.210. We have taken h/J = [panel 
(a)], h/J = 0.5 [panel (b)], h/J = 1 [panel (c)] and h/J = 5 [panel (d)]. 



C. Looking at the spin-chain as a channel 

We now slightly change perspective and consider the spins 
occupying sites j = 2, ...,N as the elements of an environment 
for spin 1 . This allows us to investigate the robustness of the 
quantum correlations shared with spin under the influence 
of an amplitude damping channel. A similar analysis can be 
found in Ref. 12711 where the entanglement dynamics in the 
presence of the same model has been studied. Here we ex- 
tend this study to the analysis of QD so as to show that, at 
variance with entanglement, a unilateral non-unitary channel 




t[J- 1 ] 



FIG. 9: (Color online) QD between spin and 3, for a chain of N=3 
spins whose pair (1,0) is initialized in an X-type state with maxi- 
mally mixed marginals having c x = 1, c y = c, = and thus 0(0) = 0. 
We have taken h/J = (blue dashed line) and h/J = 2 (full red line). 
The dot-dashed line is the transition amplitude [/i(f)|. 



can induce quantum correlations in an initial state that is fully 
classical. 

By using Eqs. (O and (O, we obtain that entanglement 
sudden death (ESD) occurs when \fi(t)\ 2 < 1 - lp ' 4 ' 2 ' P22P33 



and \Mt)\ 2 <1 



IP2_ll -P11P44 
P33P44 ' 



for p 33 p44 5t0. 



P33P44 

For states 



with initial finite entanglement, the quantities \pu\ - 
P22P33 and IP23I 2 - P11P44 are always comprised between 
and 1 . For P33 (P44) = 0, the concurrence evolves as 
C(i,0)(0 = l/i (01 C(i,0)(0). Besides this latter case, where ESD 
occurs only at \f\(t)\ = 0, mixed non-separable quantum states 
experience ESD due to the fact that appropriate choices of the 
interaction parameters make |/i (f)| range from to 1 : The am- 
plitude damping channel, acting only on one spin, is able to 
erase completely the quantum correlations which give rise to 
entanglement. 

Now let us turn our attention to the conditions for the van- 
ishing of QD discussed previously. As far as condition 1 
is concerned, we observe that a zero-QD state will remain 
such at all times because an environment addressing only 
one spin of a bipartite state cannot build up quantum coher- 
ences. If the latter are initially present in such a way to 
start from a non-zero QD state, condition 1 can be possi- 
bly fulfilled only asymptotically. On the other hand, start- 
ing from a zero-QD state according to condition 4, we note 
that the constraint p„ = 1/4 breaks down because of the dy- 
namics embodied by Eqs. © and the classical-classical state 
evolves into a quantum-classical one, so that T>~* + 0. Con- 
versely, starting from a state with non-zero QD, it is neces- 
sary, in order to satisfy condition 4, that the initial state con- 
sists of a quantum-classical state (condition 2) and that the 
environment acts on the spin that, when subjected to a com- 
plete projective measurement, gives rise to non-zero one-way 
QD. Furthermore, as the amplitude damping channel implies 
pn(0 > pn(0), only quantum-classical states with pn(0)<^ 
can evolve to a classical-classical one. This will occur at times 
f such that MO = l/2+p„/(l -2p„). 
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IV. CONCLUSIONS 

We have studied the propagation of quantum correlations 
across a chain of interacting spins by looking at the perfor- 
mance of two significant figures of merit: quantum discord 
and entanglement of formation. The amount of transported 
quantum correlations has been quantified when the chain is 
seeded with various instances of two-spin mixed states. We 
have explicitly considered the case of pure entangled states, 
as well as the members of the boundary family of MDMSs, 
which maximize QD at set values of global entropy. Discord 
appears to be consistently better transported than entangle- 
ment, except for highly pure input states whose EOF is trans- 
ported across a chain guaranteeing a large single-excitation 
transition amplitude. Moreover, we have performed a brief 
case-study on the conditions under which QD is actually cre- 
ated upon propagation across the chain, pointing towards the 
class of states giving rise to this effect, and analyzing, in par- 
ticular, how the discord increases under the effect of a mag- 
netic field. Finally we have revealed that, whereas the EoF of 
all mixed states undergoes ESD by an appropriate choice of 
the interaction parameters, QD vanishes only under very spe- 
cific initial state conditions and interaction settings. In par- 
ticular quantum-quantum states allow for vanishing discord 
under the influence of the spin environment considered here 
only asymptotically in time. 

Our work extends the investigations performed so far on 



the propagation of entanglement to the broader realm of more 
general quantum correlations proving that, in the dispersive 
medium consisting of interacting quantum particles, correla- 
tions beyond entanglement are favored. Such results motivate 
the study on a few aspects of this general problem that still 
remain to be addressed, such as the extension to other quan- 
tifier of non-classicality, their behavior under more general 
dynamical conditions of the spin media and the quantification 
of the corresponding non-Markovian effects brought about by 
the spin chain on the dynamics of general indicators of quan- 
tumness. 
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